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I.  INTRODUCTION 


The  surface  properties  of  layered  materials  such  as  graphite  and  the  chalcogenides  M0S2  and 
GaSe  are  a  subject  of  much  interest.  \V^thin  each  layer,  strong  covalent  bonding  prevails; 
however,  only  weak  van  der  Waals  interactions  are  operative  between  layers,  lb  a  first 
approximation,  the  surface  properties  of  such  materials  might  be  thought  of  as  identical  to  the 
bulk  properties.  For  instance,  the  (0001)  surfaces  of  M0S2,  NbSe2,  and  GaSe  show  1x1 
LEED  patterns,  indicative  of  no  major  surface  reconstruction.  A  careful  study  of  the  LEED 
intensities  for  M0S2  and  NbSe2  indicates  a  contraction  of  less  than  5%  of  the  first  van  der 
Waals  gap  in  M0S2,  and  no  changes  in  NbSe2  (Ref.  1). 

The  dispersion  of  surface  phonons  in  graphite  has  been  measured  by  high-resolution  electron 
energy  loss  spectroscopy  (HREELS)  ^fs.  2, 3)  and  is  considered  to  be  similar  to  that  of  the 
bulk  phonons  as  measured  by  neutron  scattering  (Ref.  4).  Inelastic  He  scattering  was  used  to 
measure  the  surface  phonon  dispersion  of  GaSe  (Ref.  5),  where  large  differences  from  the  bulk 
phonon  dispersion  were  observed.  The  lower  surface  phonon  energies  were  attributed  to 
lower  force  constants  and  higher  polarizability  of  Se  at  the  surface,  causing  increased  van  der 
Waals  attraction. 

In  this  report,  we  describe  the  surface  phonon  dispersion  of  M0S2  measured  by  HREELS. 

The  crystal  structure  of  M0S2  is  illustrated  in  Figure  la.  Each  S-Mo-S  sandwich  is  tightly 
bound  internally  and  interacts  with  neighboring  sandwiches  only  through  van  der  Waals 
forces.  Figure  lb  shows  projections  of  the  M0S2  unit  cell  on  the  (110)  plane  and  the  (0001) 
(basal)  plane,  and  the  conventional  surface  Brillouin  zone  with  the  high  symmetry  points 
marked  is  shown  in  Figure  Ic.  In  the  experiments  reported  here,  the  surface  phonons  are 
found  to  be  lower  in  energy  than  the  bulk  phonons,  as  in  the  case  of  GaSe,  indicating 
differences  between  the  bonding  in  the  topmost  layer  and  the  bulk  layers. 
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2H-MoS2 


Figure  1.  (a)  Crystal  structure  of  M0S2,  (b)  unit  cell  of  M0S2,  and  (c) 
surface  BrUlouin  zone  of  M0S2. 
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II.  EXPERIMENTAL 


Basal  plane  surfaces  of  M0S2  were  produced  by  cleavage  in  air  of  natural  molybdenite  crystals 
(supplied  by  >^^d’s  Natural  Science  Establishment,  Rochester,  NY).  The  samples  were 
cleaned  by  heating  to  750  K  for  10  min  at  10"^®  torr  to  drive  off  adsorbed  hydrocarbons.  This 
procedure  removes  virtually  all  contaminants  without  producing  defects  (Ref.  6).  HREELS 
spectra  before  cleaning  showed  only  hydrocarbon  contaminant  vibrations;  after  cleaning,  only 
those  due  to  M0S2  were  present.  Good  quality  (0001)  1x1  LEED  patterns  were  obtained  and 
used  to  orient  the  sample  along  the  <010>  or  <  110  >  azimuths  for  the  dispersion 
measurements. 

The  spectra  were  taken  with  a  double  pass  spectrometer  (LK  2000  from  LK  lechnologies,  Inc.) 
with  a  sample  current  of  - 10"^®  A  and  a  resolution  of  9  meV  for  the  M0S2  samples  used  here. 
A  resolution  of  6  meV  has  been  obtained  by  biasing  the  sample  so  that  the  electrons  are 
reflected  from  a  smooth  electric  field  above  the  sample,  but  a  real  surface  with  contaminants 
and  defects  will  have  a  poorer  resolution.  Spectra  were  collected  digitally  and  typically  took 
1.5  hr  each.  The  angle  of  incidence  of  the  electrons  was  fixed  at  60”  with  respect  to  the 
surface  normal,  and  the  angle  of  collection  of  the  scattered  electrons  was  varied  from  60” 
(specular  reflection)  to  27”  with  respect  to  the  surface  normal  to  collect  electrons  of  the 
desired  momentum  transfer.  Spectra  were  collected  with  incident  beam  energies  of  7, 15,  30, 
and  40  eV.  In  some  cases  the  sample  was  cooled  to  120  K  to  reduce  the  intensity  of  phonon 
deexitation  peaks  (“gain”  peaks). 
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III.  RESULTS 


Typical  spectra  for  15  eV  incident  beam  energy  are  shown  in  Figure  2  (along  T  to  M)  and 
Figure  3  (along  F  to  K).  Optical  phonons  near  50  meV  are  present  in  specular  reflection  and, 
as  the  momentum  transfer  is  increased,  additional  vibrations  begin  to  appear  near  the  elastic 
peak.  The  same  general  observations  can  be  made  for  the  spectra  taken  with  other  beam 
energies;  the  phonon  peaks  are  less  intense  for  the  30  and  40  eV  incident  beams.  Overtones 
of  the  vibrations  can  be  observed  at  twice  the  phonon  energy,  but  are  so  low  in  intensity  that 
they  cannot  be  used  in  the  analysis. 

The  overlap  of  the  peaks  makes  location  of  the  peak  positions  difficult,  and  a  curve-fitting 
procedure  was  used.  The  overlapping  peaks  were  represented  by  Gaussians  of  half-width 
9  meV  (the  half-width  of  the  elastic  peak  at  k  =  0).  This  half-width  was  used  for  all  the  peaks 
in  the  absence  of  any  compelling  reason  for  allowing  it  to  vary,  in  order  to  reduce  the  number 
of  adjustable  parameters  used  in  the  curve  fitting.  The  phonon  creation  and  annihilation  (loss 
and  gain)  peaks  were  equally  spaced  on  each  side  of  the  elastic  peak,  and  were  given  an 
intensity  ratio  appropriate  to  the  temperature  of  the  measurement: 

4  =  7/ 

where  Ig  and  Ii  are  the  intensities  of  the  gain  and  loss  peaks,  respectively,  E  is  the  energy  of 
the  phonon,  k  is  the  Boltzmann  constant,  and  T  is  the  temperature.  The  position  and 
amplitude  of  the  elastic  and  loss  peaks  were  then  varied  to  produce  the  best  flt  to  the  data. 

For  many  of  the  fits,  an  automated,  interactive  least-squares  procedure  was  used.  The  error 
in  peak  positions  is  estimated  to  be  ±0.2  meV  for  peaks  that  are  well  resolved  (for  instance, 
the  optical  phonons  at  around  50  to  60  meV).  However,  for  the  most  heavily  overlapped 
regions,  near  the  elastic  peak  where  five  peaks  were  sometimes  necessary  to  fit  the  data,  the 
computer  time  needed  for  this  procedure  was  excessive.  Instead,  a  variety  of  trial  fits  were 
compared  to  the  data,  and  the  best  fit  was  chosen  by  eye.  Figure  4  shows  a  typical  fit 
obtained  in  this  manner,  using  the  data  for  7  eV  incident  electrons  and  collection  at  30° 
off-specular  along  F  to  M.  From  the  comparison  process,  it  is  estimated  that  the  peak 
positions  are  reproducible  to  ±  0.5  meV  The  root-mean-square  deviation  of  the  chosen  fit 
from  the  data  was  found  to  be  less  than  that  for  the  rejected  trials  (those  with  peak  energies 
of  ±  0.5  meV  from  the  selected  values)  for  the  several  cases  where  this  deviation  was 
calculated. 

The  momentum  transfer  for  each  experiment  was  calculated  using  the  conservation  of 
momentum  parallel  to  the  surface: 

ki  =  jTmE/h  (sin0p„,  -  sin0i„c) 

where  m  is  the  electron  mass,  h  is  Planck’s  constant,  and  Oput  and  Ojnc  are  the  angles  of 
collection  and  incidence,  respectively,  measured  with  respect  to  the  surface  normal.  The  data 
are  presented  in  reduced  coordinates  (units  of  2-ir/a  for  r  to  K,  An/a  J3  for  r  to  M,  where  a  is 
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the  lattice  parameter  in  the  basal  plane)  that  reflect  the  fraction  of  the  distance  to  the  zone 
boundary.  The  dispersion  curves  are  given  in  Figure  S,  along  with  selected  calculated  bulk 
dispersion  curves  of  Wakabayashi  et  al.  (Ref.  7).  These  curves  fit  their  neutron  scattering 
measurements  of  the  bulk  phonon  dispersion  very  well,  but  exhibit  some  differences  from  the 
surface  phonon  dispersion,  as  will  be  discussed  in  section  IV. 


I _ I  V-  — -1 -  I 

-100  -50  0  50  100  150 

LOSS  ENERGY  (meV) 


Figure  2.  HREELS  spectra  of  M0S2  along  the  <  110  >  azimuth  (F  to 
K).  Incident  beam  energy  was  15  eV.  Each  spectrum  is  labeled  with 
the  angle  of  collection  with  respect  to  the  surface  normal. 
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Figure  3  .  HREELS  ^)ectra  of  M0S2  along  the  <  010  >  azimuth  (T  to 
M).  Incident  beam  energy  was  IS  eV.  Each  spectrum  is  labeled  with 
the  angle  of  collection  with  respect  to  the  surface  normal. 
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LOSS  ENERGY  (meV) 


Figure  4.  Examole  of  curve  fitting  to  HREELS  spectrum  of  M0S2  with 
7-eV  incident  beam  and  30*  collection  angle  along  <  010  >  azimuth  (f 
toM). 
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Figure  S.  Surface  phonon  dispersion  curves  for  M0S2  along  the  <  110>  and  <010>  directions.  The  lines  are  calculated 
bulk  curves  of  Wakabayashi,  et  al.  The  data  points  are  measured  for  incident  beam  energies  of  7  eV  (•),  IS  eV  (O),  40  eV 
(□,  r  to  M),  30  eV  (□,  r  to  K),  and  15  eV  using  a  cooled  sample  (A). 


IV.  DISCUSSION 


In  bulk  M0S2  there  are  six  atoms  per  unit  cell,  so  18  phonon  modes  are  expected;  3  acoustic 
modes  and  15  optical.  At  r  the  modes  can  be  labeled  by  irreducible  representations  of 
theDSi  point  group  (P3  /mmc  space  group).  These  modes  were  listed  by  Agnihotri  et  al.  (Ref. 
8),  along  with  diagrams  of  atomic  dispacements  for  each  mode.  The  unit  cell  of  the  surface 
layer  of  M0S2  has  three  atoms  (two  sulfur  and  one  molybdenum),  so  9  modes  are  expected:  3 
acoustic  -nodes  and  6  optical  modes.  The  point  group  for  the  surface  is  Ccv;  there  will  be  3 
Al  modes  and  3  doubly  degenerate  E  modes.  These  modes  are  listed  in  Table  I,  with  their 
related  bulk  modes  (Ref.  8),  atomic  displacements  (following  Agnihotri  et  al.),  related 
measiured  bulk  phonon  energies,  and  the  surface  phonon  energies  at  F  determined  in  this 
work. 


llible  I.  Lattice  Modes  at  F  for  M0S2 


IRREDUCIBLE 

REPRESENTATION 

Ai 
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A, 
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Al 
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RELATED  BULK 
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REPRESENTATIONS 
(ref  8) 

^19 

^2u 

aL 

cl 

‘'2g 

^2u 

Eig 

ATOMIC 

DISPLACEMENTS 

6 

Q* 

• 

(i 

•*Q 

(2K 

Q* 

BULK  ENERGY 
(meV)  (ref  7) 

50.7 

35.6 

57.8 

46  7 

0 

0 

SURFACE  ENERGY 
(meV) 

47.6 

57.2 

- 

0 

0 

^PROJECTED ON (110) PLANE.  Q=S,  •=Mo,  <0001> 
DIRECTION  IS  TOWARD  THE  TOP  OF  THE  PAGE 


The  HREELS  peaks  measured  in  these  experiments  can  be  assigned  to  the  expected  surface 
phonon  modes.  Two  of  the  Ai  modes  are  optical  and  will  be  dipole-active  vibrations,  lliese 
should  be  observed  in  specular  reflection,  and  should  decrease  in  intensity  as  the  elastic  peak 
decreases  in  intensity  when  scattering  is  detected  off-specular.  The  peaks  at  47.6  and 
57.2  meV  behave  in  this  manner,  so  they  are  assigned  as  Ai. 
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There  are  two  E  optical  modes,  which  are  in-plane  vibrations  and  could  be  excited  by  the 
impact  scattering  mechanism  that  dominates  HREELS  experiments  away  from  the  specular 
direction.  For  our  experimental  geometry,  the  observed  impact  modes  must  be  totally 
symmetric  with  respect  to  the  scattering  plane.  These  two  modes  would  then  be  visible  in  the 
<  110  >  direction,  and  along  this  azimuth  we  have  a  peak  in  the  spectra  for  k  >  0.2.  The  E 
modes  should  not  be  visible  along  the  <010>  direction,  and  as  expected  we  find  no 
corresponding  peak  in  the  spectra.  The  energy  of  the  peak  found  along  the  <  110>  direction 
is  not  listed  in  Table  I  because  it  was  not  measured  at  F. 

Acoustic  modes  approach  0  meV  as  k  approaches  0.  We  expect  an  Ai  mode  (the  Rayleigh 
wave)  to  be  the  lowest  energy  vibration,  and  assign  the  lowest  energy  peak  in  the  spectra  to 
this  mode.  The  acoustic  E  mode  is  present  only  along  the  <010>  azimuth. 

At  k  =  0.5  to  0.6  along  the  <  110  >  azimuth,  two  peaks  with  very  low  signal  to  noise  are 
found.  The  assignment  of  these  peaks  is  unclear. 

The  energies  of  the  surface  phonons  are  less  than  those  of  the  related  bulk  phonons.  The  Aj 
optical  modes  are  0.6  and  3.1  meV  lower  in  energy  (experimental  error  ±0.2  meV),  with  a 
similar  dispersion  to  the  bulk  modes.  The  two  E  optical  modes  are  probably  degenerate,  since 
we  only  observe  one  peak  in  the  spectra  that  we  can  assign  to  modes  of  this  type.  They  must 
also  be  lower  in  energy  than  their  related  bulk  modes,  which  were  not  measured  along  F  to  K, 
since  Rayleigh’s  theorem  states  that  if  the  spring  constant  of  one  vibration  in  a  set  of 
harmonically  coupled  masses  is  lowered,  all  the  vibrational  frequencies  must  decrease.  The 
energies  of  the  acoustic  modes  measured  here  are  generally  30%-50%  lower  than  those  of  the 
bulk  acoustic  modes  measured  by  Wakabayashi  et  al.  (Ref.  7)  across  the  Brillouin  zone.  These 
differences  are  larger  than  our  estimated  experimental  error  ( ±  0.5  meV). 

The  large  softening  of  the  surface  phonons  with  respect  to  the  bulk  phonons  is  similar  to  the 
case  for  GaSe,  and  different  from  that  of  graphite.  The  optical  phonons  in  GaSe  were  also 
softened  by  several  percent.  The  explanation  advanced  for  GaSe  surface  phonon  softening 
could  also  be  viable  for  M0S2:  increased  S  polarizability  at  the  surface,  and  increased  van  der 
Waals  attraction  at  the  first  gap.  This  explanation  is  also  supported  by  the  LEED  results, 
indicating  a  contraction  of  the  first  van  der  Waals  gap  in  M0S2. 
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V.  CONCLUSIONS 


The  dispersion  of  the  surface  phonons  of  M0S2  has  been  measured  by  HREELS  along  both 
the  <010>  and  <  110  >  azimuths  (T  to  K  and  F  to  M).  The  surface  phonons  are  lower  in 
energy  than  the  corresponding  bulk  phonons,  as  was  found  in  previous  studies  of  GaSe.  The 
contraction  of  the  first  van  der  Waals  gap  in  M0S2  as  measured  by  LEED,  combined  with  the 
softer  surface  phonons  measured  in  these  experiments,  suggests  that  the  S  atoms  at  the 
surface  may  be  more  polarizable  than  those  in  the  bulk,  as  seems  to  be  the  case  for  the 
surface  Se  atoms  in  GaSe. 
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